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ABSTRACT 



The magnetic field in stellar radiation zones can play an important role in phenomena such as mixing, angular momentum transport, 
etc. We study the elfect of rotation on the stability of a predominantly toroidal magnetic field in the radiation zone. In particular we 
considered the stability in spherical geometry by means of a linear analysis in the Boussinesq approximation. It is found that the 
effect of rotation on the stability depends on a magnetic configuration. If the toroidal field increases with the spherical radius, the 
instability cannot be suppressed entirely even by a very fast rotation. Rotation can only decrease the growth rate of instability. If the 
field decreases with the radius, the instability has a threshold and can be completey suppressed. 
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1. Introduction 

It is rather uncertain which magnetic field can be present in 
stellar radiation zones but this field can play an important role 
in many phenomena in stars such as mixing, angular momen- 
tum transport, formation of tachocline etc. (see, e.g., Gough 
& Mclntyre 1998, Heger et al. 2005, Eggenberger et al. 2005, 
Mathis & Zahn 2005). Likely, dynamos cannot operate in radia- 
tion zones, where no strong flows are available to sustain a vig- 
orous dynamo action. Possibly, relic magnetic fields acquired by 
the star at the early stage of evolution can persist there. This type 
of fields could have formed, for instance, because of differential 
rotation that could have stretched the lines of a weak primordial 
seed field into a dominant toroidal field. Since the conductivity 
is high, the large-scale relic field could survive in the radiation 
zone during the life-time of a star. 

The magnetic field, however, can evolve in a radiation zone 
also because of the development of various instabilities. For ex- 
ample, the magnetorotational instability could occur if the radia- 
tion zone is magnetized and rotates differentially. However, dif- 
ferential rotation is unlikely in radiation zones and, perhaps, can 
exist only in stellar tachoclines. Over the past decade, instabili- 
ties of the stellar tachocline have been extensively studied (see, 
e.g., Dikpati et al.(2009) and reference therein). The tachocline 
is thin and its stability properties are rather peculiar. For in- 
stance. Oilman & Fox (1997) showed that the tachocline lati- 
tudinal shear is unstable to nonaxisymmetric disturbances when 
a toroidal magnetic field is present. Instabilities in the tachocline 
have been studied in detail by Dikpati et al.(2009) for a wide 
range of rotation and toroidal field profiles. Since rotation is rigid 
in radiation zones, instabilities of the magnetic field most likely 
are cuiTent-driven. These instabilities do not require differential 
rotation, and they are well studied in cylindrical geometry in the 
context of laboratory fusion research (see, e.g., Freidberg 1973, 



Ooedboed 1971, Ooedbloed & Hagebeuk 1972, Tayler 1973a,b, 
1980). In astrophysical conditions, the instability caused by elec- 
tric currents might have a number of characteristic features even 
in cylindrical geometry (see Bonanno & Urpin 2008a,b, 2011). 
The nonlinear evolution of the Tayler instability was studied by 
Bonanno et al. (2012), who argued that symmetry-breaking can 
give lise to a saturated state with nonzero helicity even if the ini- 
tial state has zero helicity. A production of nonzero heficity is 
crucial for the dynamo action. 

The stability of the spherical magnetic configurations is stud- 
ied in much less detail. This problem is of particular interest in 
relation to Ap star magnetism (Braithwaite & Spruit 2004). With 
numerical simulations Braithwaite & Nordlund (2006) studied 
the stability of a random initial field. They argued that this 
field relaxes on a stable mixed magnetic configuration with both 
poloidal and toroidal components. A study of the magnetic con- 
figurations with a predominantly toroidal field is of particular 
importance for radiation zones because this field can be easily 
formed by differential rotation at the early evolutionary stage. 
Numerical modeling by Braithwaite (2006) confirmed that the 
toroidal field with oc s or oc (s is the cylindrical radius) is 
unstable to the m = 1 mode, as was predicted by Tayler (1973a). 
The stability of azimuthal fields has also been studied by Spruit 
(1999). The author used a heuristic approach to estimate the 
growth rate and criteria of instability. Unfortunately, many of 
these estimates and criteria are valid only near the rotation axis 
and do not apply in the main fraction of the volume of a radiation 
zone where the stability properties can be qualitatively different 
(see Zahn et al. 2007). Recently, Bonanno & Urpin (2012) have 
considered the stability of the toroidal field in radiation zones by 
making use of a linear analysis and taking into account stratifica- 
tion and thermal conductivity. It is widely believed that stratifica- 
tion can suppress the Tayler instability. Bonanno & Urpin (2012) 
calculated the growth rate of instability and argued that the stabi- 
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lizing influence of gravity can never entirely suppress the insta- 
bility caused by electric currents. However, a stable stratification 
can essentially decrease the growth rate of instability. 

In this paper, we consider the effect of rotation on the sta- 
bility of magnetic configurations with a predominantly toroidal 
field. Rotation is often considered as one more factor that can 
suppress the Tayler instability and stabilize the magnetic con- 
figurations. For instance. Spruit (1999) found that the growth 
rate of the Tayler instability in a rotating star should be on the 
order of ~ u)A{ii>Al^) if O » u)a, where oja and O are the 
Alfven frequency and angular velocity of the star, respectively. 
StabiUty of the toroidal field in rotating stars has been consid- 
ered by Kitchatinov (2008), and Kitchatinov & Riidiger (2008), 
who argued that the magnetic instability is determined by the 
threshold field strength at which the instabihty sets. Estimating 
this threshold in the solar radiation zone, the authors imposed the 
upper limit on the magnetic field to be ^ 600 G. The stability of 
the toroidal field in a rotating radiation zone has been studied by 
Zahn et al. (2007) in the particular case oc s. The particular 
type of oscillatory modes found by these authors is related to 
rotation and is stable in the nondissipative Umit. However, insta- 
bility can occur in the form of an oscillatory diffusive instability 
if dissipation is provided by radiative or Ohmic diffusion. 

The paper is organized as follows. The basic equations and 
mathematical formulation of the problem are presented in Sec. 2. 
This is followed by results of numerical calculations of the 
growth rate and frequency of the instability in Sec. 3. The pa- 
per closes with a summary of the main results and some remarks 
in Sec.4. 

2. Basic equations 

Consider the stabihty of an axisymmetric toroidal magnetic field 
in the radiation zone using a high conductivity limit. We work in 
spherical coordinates (r, 6, (p) with the unit vectors (e^, eg, e^). 
We assume that the radiation zone rotates with the angular ve- 
locity Q=const and that the toroidal field depends on r and 6, 
= B^{r, 6). If the magnetic field is subthermal (so that the 
magnetic pressure is lower than the gas pressure), one can apply 
the incompressible limit for a consideration of low-frequency 
modes. We conduct the analysis of hydromagnetic modes in the 
rest frame rather than in the corotating frame. GeneraUy, some 
rotation-related effects can be missing in this case. Some in- 
sights into how hydromagnetic waves behave in the corotating 
frame can be found in the papers by Hide (1969) and Acheson 
& Hide (1973). For instance, effects due to the density strati- 
fication can be neglected when considering the dispersion re- 
lation for hydromagnetic osciUations if the Brunt- Vaisala fre- 
quency is much lower than twice the angular velocity. However, 
the Brunt- Vaisala frequency is basically much higher than the 
angular velocity in stellar radiation zones and these effects can 
be neglected. 

The MHD equations in the incompressible Umit read 



-I- (v • V)v 
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where g is gravity. The equation of thermal balance reads in the 
Boussinesq approximation 



where V^^r is the adiabatic temperature gradient. 

In the basic (unperturbed) state, the gas is assumed to be in 
hydrostatic equilibrium, then 



— = g+—iVxB)xB+esQ^ rsinO, 
p Anp 



(5) 



dT 

-.v.(vr. 



v„dr) = 0, 



(4) 



where e , is the unit vector in the cylindrical radial direction. The 
rotational energy is assumed to be much lower than the gravita- 
tional one, g » rQ^. The origin and structure of the magnetic 
field in radiation zones are unknown. However, g is approxi- 
mately radial in our basic state since we assume that the mag- 
netic energy is subthermal. Only small variations of the density 
and temperature are required in the meridional direction to bal- 
ance the centrifugal and Lorentz forces for a given magnetic con- 
figuration. 

We consider a linear stability. Weak perturbations will be in- 
dicated by subscript 1, while unperturbed quantities will have 
no subscript. Linearizing Eqs.(l)-(4), we take into account 
that weak perturbations of the density and temperature in the 
Boussinesq approximation are related by pi/p - -/3{Ti/T), 
where /3 is the thermal expansion coefficient. For weak perturba- 
tions, we use a local approximation in the 0-direction and assume 
that their dependence on 6 is proportional to exp(-;70), where 
/ » 1 is the longitudinal wavenumber. Since the basic state is 
stationary and axisymmetric, the dependence of perturbations on 
t and ip can be taken in the exponential form as weU. Then, per- 
turbations are proportional to exp (erf - ilO - imcp), where m is 
the azimuthal wavenumber. The corresponding wavevectors are 
kg = l/r and = m/r sin 6, respectively. The dependence on r 
should be determined from Eqs.(l)-(4). 

The problem of the magnetic field stability is complicated 
from both the physical and computational points of view. The 
conclusions of different studies are often contradictory (com- 
pare, for example, the results by Spruit (1999) and Kitchatinov 
& Ruediger (2008)), and very often the authors do not discuss in 
detail the reason of the controversy. Because of the complexity 
of this problem, in our opinion, the best approach is to separately 
consider the role played by different physical factors (gravity, ro- 
tation, conductivity, etc.) in modifying the global structure of the 
unstable modes. Only after this process has been clarified is it 
possible to arrive at a global picture of the parameter space. The 
main motivation of this paper is therefore to clarify the effect 
of rotation by considering different radial profiles of the toroidal 
field, in contrast to the study proposed in detail by Bonanno & 
Urpin (2012), where only the effect of gravity was considered. 
Therefore, we consider a simplified problem assuming that strat- 
ification is neutral and Vr = ^adT- In this case, the stabilizing 
effect of gravity is neglected and we can study how the instabihty 
is affected by rotation alone. For the sake of simplicity, we also 
assume that the unperturbed density is approximately homoge- 
neous in the radiation zone. As we mentioned above, meridional 
variations of the density should be small in a subthermal mag- 
netic field. A radial variation is not small in real stars but omit- 
ting it does not qualitatively change the main conclusions and 
substantially simplifies calculations. Eliminating all variables in 
favor of vir, we obtain with the accuracy in terms of the lowest 
order in ikgr)~^ 

(cr^ + 0)1 + £)Qf) v;; + 1^0-2 + ^a^lj (6) 

^crl-kl(crl + col)-DQ.lkl+ 
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-CO. 



2kl 



D - i(TQ^e 



r rk^ (7^ 



Vlr = 0, 



where the prime denotes a derivative with respect to r and 



(Tq — cr — imO., 



2 '^v"f 

47rp 



D = 



(7) 



i2i = 2Qcos6l, Q<, = 2Qsin6l, kl ^ kj + kl, 

^ = |lnW. 

If Q = 0, Eq.(6) transforms into the equation derived by 
Bonanno & Urpin (2012). 

Some important stability properties of the toroidal field can 
be derived directly from Eq.(6). Consider perturbations with a 
very short radial wavelength for which one can use a local ap- 
proximation in the radial direction, such as v\r °^ exp{-ikrr), 
where kr is the radial wavevector. If A:, » max(ke,k^), then 
Eq.(6) yields with the accuracy in terms of the lowest order in 
(ferO"' the following dispersion equation 



or 



0, 



(8) 



(9) 



It is easy to show that this equation has only imaginary roots. 
Therefore, modes with a short radial wavelength are always sta- 
ble to the current-driven instability in contrast to the result ob- 
tained by Kichatinov (2008) and Kichatinov & Rudiger (2008). 

3. Numerical results 

We assume that the radiation zone is located at Ri < r < R or, 
introducing the dimensionless radius x = r/R, at x, < x < 1 
where x, - R JR. We choose the intemal radius of the radiation 
zone, Xi, to be equal to 0. 1 from computational reasons. We have 
verified that our results are basically insensitive to the precise 
value of X, as long as it is close to the center. 
The toroidal field can be represented as 



- Bi)^{x) sin 9, 



(10) 



where Bq is the characteristic field strength and ^ ~ 1 is a func- 
tion of the spherical radius alone. The dependence of i/f on x 
is uncertain in the radiation zone and, in this paper, we con- 
sider three different possibilities. Fig. 1 plots the profiles if/{x) 
for the models (1), (2), and (3) used in our calculations. The 
case where the field reaches its maximum at the outer bound- 
ary (model (2)) can mimic, for example, the radiation zone of 
a star with a convective envelope. In this case, the bottom of a 
convection zone likely is the location of the toroidal field gen- 
erated by a dynamo action. The toroidal field can penetrate into 
the radiation zone, for instance, because of diffusion. Model (2) 
can also mimic the toroidal field in the liquid core of neutron 
stars. Likely, the magnetic field of these objects is generated by 
turbulent dynamo during the very early phase of evolution when 
the neutron star is subject to hydrodynamical instabilities (see 
Bonanno et al. 2005, 2006). A large-scale dynamo is most efli- 
cient in the surface layers where the density gradient is highest. 
Therefore, the generated field increases outward and reaches its 
maximum in the outer layers (Bonanno et al. 2005, 2006). This 
magnetic field can be subject to current-driven instabilities af- 
ter the end of the initial phase. Note that dynamo induced by 




Fig. 1. Dependence of the toroidal field on the spherical radius 
for models (1), (2), and (3). 



turbulent motions generates not only a large-scale field but also 
small-scale fields of complex topology (Urpin &, Gil 2004). The 
stability properties of the configurations with mixed small- and 
large-scale fields can be very particular and are not studied yet. 
The model (3) with a decreasing toroidal field can mimic a star 
whose magnetic field is generated in the inner convective core. 
Generally, all three models (1), (2), and (3) can be representative 
of the stars with relic magnetic fields because the details of the 
formation of these fields are very uncertain. 
Introducing the dimensionless quantities as 



r = 



T] = 



2Q 



(11) 



where aJ^^^ = B^/AnpR^, we can transform Eq.(6) into a dimen- 
sionless form. This equation with the corresponding boundary 
conditions describes the stability problem as a nonlinear eigen- 
value problem. Fortunately, the main quaUtative features of this 
problem are not sensitive to the choice of boundary conditions. 
That is why we choose the simplest boundary conditions and as- 
sume that Vlr = at r = Ri and r = R. Generally, solutions of 
Eq.(6) are complex. It is more convenient to split all quantities 
into the real and imaginary parts and to solve numerically the 
set of two real coupled equations that follows from Eq.(6). Note 
that coefficients of Eq.(6) and the corresponding dimensionless 
equation depend on 0, which in turn leads to the dependence of 
r on it. 

The stability properties in the spherical geometry are qualite- 
tively different from those in the cylindrical geometry (see 
Bonanno & Urpin 2012). Therefore, the results obtained, for in- 
stance, for the toroidal field dependending on the cyhndrical ra- 
dius alone (see Spruit 1999, Zhang et al. 2003) does not apply 
to more general magnetic configurations. The stability problem 
becomes particularly complex if the radiation zone rotates . 

In Fig. 2, we plot the growth rate and frequency of the Tayler 
modes as functions of 77 at different 6 for the model (2). Like the 
case of a non-rotating star, the Tayler instability is the most effi- 
cient at the equator (see also Bonanno & Urpin 2012). At small 
77, the growth rate is of the order of 1 Acoao but it clearly shows 
some suppression for a faster rotation. Suppression becomes sig- 
nificant already at relatively low values of 77 ~ 2 - 3. The growth 
rate decreases with an increase of 77 approximately as 1 Irj and 
it does not vanish even at very large 77. A similar behavior was 
obtained by Spruit (1999), who considered the instabihty near 
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Fig. 2. Growth rate (left panel) and frequency (right panel) of 
the Tayler modes as functions of the rotational parameter r/ for 
6 = 7t/2 (soUd), 7t/6 (dashed), and n/W (dash-and-dotted) and 
for model (1). The longitudinal and azimuthal wavenumbers are 
/ = 10 and m= I, respectively. 



the rotation axis for - B^(s). It turns out that rotation can 
never entirely suppress the Tayler instability of model (2) but 
only decreases the growth rate. Tayler modes are oscillatory in 
a rotating radiation zone in contrast to the nonrotating case. The 
frequency is basically comparable to the growth rate and also 
decreases when the rotation becomes faster. 

Fig. 3 shows the results of calculations of the growth rate and 
frequency for model (1). The results are very similar to those 
obtained for model (2). As in the previous case, the instability is 
most efficient at the equator and its growth rate decreases when 
it approaches the rotation axis. This result is at variance with 
the widely accepted opinion that toroidal magnetic configura- 
tions are always unstable at the axis (see, e.g.. Spruit 1999). This 
opinion is usually based on the similarity of the spherical mag- 
netic configuration near the axis and the axisymmetric cylindri- 
cal configuration. However, this analogy is generally incorrect 
because in spherical geometry, the toroidal field near the axis 
also depends on the radial coordinate along the axis. Therefore, 
it can be unjustified to apply the results obtained for a cyhn- 
der with - B^(s) to the case of plasma near the symmetry 
axis in spherical geometry (Bonanno & Urpin 2012). Indeed, 
in the latter case stability can crucially depend on the profile 
of the toroidal field along the symmetry axis. Like model (2), 
the instability of model (1) is strongly suppressed by rotation. 
Suppression becomes essential at t] > 3 - 4. At large rj, the 
growth rate decreases oc I/;;, as in the previous case. The sta- 
bilization cannot be reached at any large 77 but the growth rate 
can be substantially reduced. 

In Fig. 4, we plot the growth rate and frequency as func- 
tions of 77 for model 3. Magnetic configurations with a rapidly 
decreasing toroidal field are stable in a cylindrical geometry. For 
example, stability properties of the toroidal field are determined 
by the parameter a = dlnB^(s)/dlns. The field is unstable to 
axisymmetric perturbations if a > 1 and to nonaxisymmetric 
perturbations if a > -1/2 (Tayler 1973a,b, 1980). The situ- 
ation is qualitatively different in spherical geometry and even 
the toroidal field, decreasing rapidly with r, can be unstable. 
Again, the instability is most efficient at the equator. However, 
its efficiency decreases drastically toward the pole and the in- 
stability does not occur in the region around the rotation axis. 



Fig. 3. Same as in Fig.2 but for model (2). 




Fig. 4. Same as in Fig.2 but for model (3) and - n/2 (sohd) 
and ;r/3 (dashed). 

The critical angle, 6cr, which distinguishes the stable and unsta- 
ble regions, is ~ 40°. At variance with the models (1) and (2), 
the instability of the model (3) is determined by the threshold 
field strength. The threshold is not very high and corresponds to 
77 « 7. Therefore, the Tayler instability is entirely suppressed in 
model (3) if f2 > 3.5wao- The growth rate is vanishing every- 
where in the radiation zone for more rapidly rotating stars. 

4. Conclusion 

We have considered the stability of the toroidal field in rotat- 
ing stellar radiation zones. The stabihty properties of the spheri- 
cal magnetic configurations turn out to be qualitatively difl'erent 
from those of the cylindrical configurations (Bonanno & Urpin 
2012). For instance, in contrast to widely accepted opinion, the 
toroidal field can be stable even near the symmetry axis. The 
reason for this behavior is the spatial dependence of the mag- 
netic field along the axis which can provide a stabilizing efl'ect. 
Therefore, a direct analogy between the stability of a cylinder 
with the azimuthal field and the toroidal field near the axis in 
radiation zones is generally incorrect. 

Rotation provides a stabilizing eflfect on the Tayler instabil- 
ity. The effect of rotation can be characterized by the parameter 
T], which is generally large in radiation zones. It turns out that 
the effect of rotation depends critically on the magnetic config- 
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uration. If the toroidal field increases with the spherical radius 
within the radiation zone (or some fraction of it, see models (1) 
and (2)) then rotation cannot entirely suppress the Tayler insta- 
bility even at very large 77. The instability growth rate is dis- 
cernible for any rotation but can be substantially reduced at very 
large 77. A reduction of the growth rate becomes important even 
at a relatively low angular velocity Q ~ 2 - 3a>Ao- At large 77, 
the growth rate behaves like cijao((^ao/^), as was obtained by 
Spruit (1999) for a particular case = B^p(s) near the rotation 
axis. The Alfven timescale, cj^q, is short compared to the stellar 
life-time, therefore even a suppressed instability with a reduced 
growth rate can be significant for radiation zones. It should be 
noted also that, most likely, the field does not decay to zero 
because of this instability. When the field becomes weaker, the 
growth rate of the instability decreases and the field cannot de- 
cay to values lower than those resulting from the condition that 
the growth rate is on the order of the inverse life-time of a star. 
Therefore, a weak field can change only insignificantly during 
the Ufe of the star, although its radial profile can be unstable. 

If the toroidal field decreases with the spherical radius (see 
Fig. 4 for model (3)), the rotation efl'ect is different. In contrast to 
models (1) and (2), the instability of model (3) is determined by 
the threshold field strength. The threshold is not very high and 
corresponds to 77 w 7. Therefore, the Tayler instability is entirely 
suppressed in the model (3) if£l> S.Scoao and modes are stable 
everywhere in the radiation zone for more rapidly rotating stars. 
Higher eigenmodes are suppressed more strongly than the fun- 
damental one and perturbations with a short radial wavelength 
are always stable. Since instability is not suppressed atrj < 1, 
this implies that the magnetic field should satisfy the condition 
Bq > ^lAnp. Estimating 0/^-2x10^ cm/s and p ~ 0.1 

g/cm^, we obtain that instability can arise in the radiation zone of 
the Sun if Sq ^ 7 X 10"* G. This estimate is more than two orders 
of magnitude higher than that obtained by Kichatinov & Riidiger 
(2008). These authors considered stability of the toroidal field 
assuming that perturbations are global in the meridional direc- 
tion and short-scaled in radius. As a result, they obtained that the 
most rapidly growing modes modes are indefinitely short in the 
radial direction if diffusion is neglected. This conclusion is ques- 
tionable because simple analytic considerations (see Section 2 
of the present paper) show that perturbations with short radial 
wavelength should always be stable. In particular our numeri- 
cal calculations also show that the fundamental eigenmode has a 
higher growth rate than higher eigenmodes. 

Finally, we calculated the growth rate assuming a neutral 
stratification of the radiation zone. If the stratification is sta- 
ble, gravity provides an additional stabilizing influence on the 
Tayler instability, as was argued by Bonanno & Urpin (2012). 
Therefore, suppression of the instability can be even more pro- 
nounced and the calculated growth rates are likely to be only 
upper Umits in real, magnetized stellar interiors. 
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